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Abstract 

The density of states of the two-dimensional fermionic Hubbard model 
in the perpendicular homogeneous magnetic field is calculated using the 
strong coupling diagram technique. The density of states at the Fermi 
level as a function of the inverse magnetic induction oscillates, and the 
frequency of these oscillations increases by an order of magnitude with the 
change of the deviation from half-filling from small to moderate values. 
This frequency variation is caused by the change of Landau subbands 
contributing to the density - in the former case they are at the periphery of 
the Landau spectrum, while in the latter case the dominant contribution is 
provided by bands near its center. With changing induction these groups 
of bands behave differently. For small deviations from half-filling the 
calculated oscillation frequency is comparable to that observed in quantum 
oscillation experiments in yttrium cuprates. 


1 Introduction 

Theoretical investigations of systems of strongly correlated electrons in strong 
magnetic fields were started shortly after the discovery of the high-Tc super¬ 
conductivity. A number of works was carried out on small clusters using the 
exact diagonalization (see, e.g., Refs. dll IS)- It is worth noting that due to 
the Peierls factor [1] , the translation symmetry of the system is changed [S] - in 
moderate magnetic fields the size of the elementary cell increases significantly. 
Clusters with sizes smaller than the size of this supercell violate the symmetry 
of the Hamiltonian and, therefore, it is difficult to extend the obtained results 
to larger crystals. Another approach used for this problem is the mean-field 
approximation (see, e.g., Refs. [S1I71IE])- The main shortcoming of this approx¬ 
imation is the neglect of the dynamic character of strong correlations. 

The interest in this problem was revived with the observation of low-fre¬ 
quency quantum oscillations in the mixed state of underdoped yttrium cuprates 
[s uni EH [ig. Based on the Onsager-Lifshitz-Kosevich theory for metals [13] 
the decreased quantum oscillation frequencies were interpreted as a manifesta¬ 
tion of small Fermi surface pockets M- To explain the appearance of these 
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small pockets proposals for various states with broken translational symmetry 
were suggested [niiiniin]. Other theories for explaining the decreased quan¬ 
tum oscillation frequency suppose that it is connected with superconducting 
fluctuations [mils] or use phenomenology of the marginal Fermi liquid m- 
Crystals, in which the decreased quantum oscillation frequencies were ob¬ 
served, are underdoped cuprates, and, therefore, they are characterized by 
strong electron correlations. Theoretically the behavior of such crystals in strong 
magnetic fields is poorly known. In this work we use the strong coupling diagram 
technique EH HI 12311^113 US] for investigating the density of states (DOS) of 
the two-dimensional (2D) fermionic Hubbard model in a perpendicular homoge¬ 
neous magnetic field. This approach allows us to consider large enough crystals 
in fields of moderate intensities. As known [illTlIlHj, the energy spectrum of 
a weakly correlated metal consists of the Landau subbands, which appear in 
the crystal potential instead of the Landau levels of free electrons. We found 
that in the approximation of a local self-energy/irreducible part each Landau 
subband forms strongly correlated bands independently of other subbands. Us¬ 
ing for the irreducible part the Hubbard-I approximation [29] we revealed that 
the DOS at the Fermi level oscillates with changing inverse magnetic induction. 
The frequency F of these oscillations increases by an order of magnitude when 
the deviation of the electron filling h from half-filling grows from a few percent 
to 25%. For the Hubbard repulsion U = 8t, t being the hopping constant, and 
the intersite distance a = dA the obtained frequency is of the order of 1 kT for 
small values of |1 — fi|. This frequency is close to those observed in the men¬ 
tioned experiments on quantum oscillations. The origin of the strong variation 
of the frequency with the electron filling is in the difference in Landau subbands 
contributing to the DOS at the Fermi level. In the case of small deviations from 
half-filling these subbands are located at the periphery of the spectrum, while 
for larger values of |1 — h| the subbands near the central part of the spectrum 
make the main contribution. The behavior of these two groups of subbands with 
the change of the magnetic induction is different. The influence of the Zeeman 
term on the DOS oscillations are also considered. 


2 Main formulas 


The Hamiltonian of the Hubbard model in the magnetic field reads 


H 


^ hv exp I A(r)dr^ a\,^aia 
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where 2D vectors 1 and 1 ' label sites of a square plane lattice, tr = ±I is the 
projection of the hole spin, and ai^ are electron creation and annihilation 
operators. The first, kinetic, term of the Hamiltonian contains the hopping ma¬ 
trix element tip and the exponential factor with the Peierls phase [3], in which 
A(r) is the vector potential. The second, Zeeman, term of the Hamiltonian 
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contains the g-factor g Ki 2, the Bohr magneton /re and the magnetic induc¬ 
tion B of the external magnetic field. It is supposed that the field is directed 
perpendicularly to the model plane, homogeneous and is only weakly disturbed 
by internal currents m- The last term of Hamiltonian m combines the on¬ 
site Coulomb repulsion with the Hubbard constant U and the electron number 
operator nic, = al^aio- as well as the term with the chemical potential fi. 

In the following consideration the Landau gauge is used, in which A(l) = 
—BlyX, where ly is the y component of the site vector 1 and x is the unit vector 
along the x axis. If we suppose that only the nearest neighbor hopping constant 
is nonzero, tu' = t "''^here a are four vectors connecting nearest neigh¬ 

bor sites, the Peierls exponential in the kinetic term of the Hamiltonian can be 
written as 

e*'^-', Ka = (2) 

n 

where ax is the x component of the vector a and y is the unit vector along the 
y axis. 

Our consideration will be restricted to the fields satisfying the condition 

e n n' 

—Ba^ = 27r—, (3) 

h n 

where a = |a|, n and n! < n are integers with no common factor. In this 
case the kinetic term of Hamiltonian © defines its translation properties - the 
Hamiltonian is invariant with respect to translations by the lattice period along 
the X axis and by n lattice periods along the y axis. To retain this symmetry 
we apply the periodic Born-von Karman boundary conditions to the sample 
with Nx sites along the x axis and nNy sites along the y axis. The boundary 

conditions define the set of allowed wave vectors with integer Ux 

and ny. As can be seen from m and o, the momenta /«a belong to this set of 
allowed wave vectors. 

Let us consider the electron Green’s function 

G'^(lV,lr) = {raiv(r)ai.(r)), (4) 

where the statistical averaging denoted by the angular brackets and time de¬ 
pendencies of the operators 

aicr(T) = and aicr(r) = 

are determined by Hamiltonian o, B is the time-ordering operator which ar¬ 
ranges operators from right to left in ascending order of times r. Hamiltonian 
o retains the spin projection and, therefore, Green’s function (j?]) is diagonal in 
this parameter. To calculate this function we use the strong coupling diagram 
technique [n [12 Ha [22 Ha [26], in which it is presented as the serial expansion 
with the unperturbed Hamiltonian Hq given by the last term of Hamiltonian 
©. In the considered case the role of perturbation, over which the power ex¬ 
pansion is carried out, is played by the first two terms of this Hamiltonian. For 
brevity, the sum of this two terms is denoted as 

ll'cr 
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Terms of the serial expansion are constructed from on-site cumulants of the 
electron operators a\a-{T) and aia-{T) and hopping lines corresponding to the 
Hamiltonian Hi (though the Zeeman term in Hi does not lead to the transfer 
of an electron to another site, we retain the term ’’hopping line“ used in this 
diagram technique). The averaging and time dependencies of operators in the 
cumulants are determined by the Hamiltonian Hq . As in the diagram technique 
with the expansion in powers of an interaction, in the strong coupling diagram 
technique the linked-cluster theorem allows one to discard disconnected dia¬ 
grams and to carry out partial summations in remaining connected diagrams. 

The diagram is said to be an irreducible one if it cannot be divided into two 
disconnected parts by cutting some hopping line. The sum of all irreducible 
diagrams without external ends is termed the irreducible part Ar°’(l'T', Ir). In 
terms of this quantity the equation for Green’s function reads 

G"(lV,lr) = if'^(lV,lr)+^ /V-(lV,liTi) 

hi; -^0 

xT‘^(lill)G‘^(i;Ti,lr)dri, (5) 


or after the Fourier transformations over space and time variables 

G'^(q'qm) = Ar'^(q'qm) -|- ^ Ar'^(q'q;^TO) 

qiq'i 

xT'^(qiq'i)G‘"(q'iqm), (6) 


where (3 = 1/T is the inverse temperature, m is an integer determining the 
Matsubara frequency LOm = (2m -I- l)7rT and 


T"(qq') = t e*«><5q,q- + ,„.+e-*«>dq,q- 


-1-2 COs((3'ya)(5qq/ -f -g/isHcr^qq/, 


(7) 


Qx and Qy are components of the wave vector q, which belong to the first Bril- 
louin zone, and k, = ^y. In the derivation of (O we took into account that /«a 
in m belongs to the set of momenta determined by the chosen periodic bound¬ 
ary conditions. Equations ([S|) and are forms of the Larkin equation m- 
Analogous equations were used for calculating the electron Green’s function of 
the Hubbard model in the absence of the magnetic field [nj mm [Ml nulls]- 
In (O, G‘^(q'qm) and A'°’(q'qm) are not diagonal with respect to momenta 
due to the reduced translation symmetry of the problem. However, from symme¬ 
try arguments one can see that these quantities are nonzero only for q^ = q -I- 
with an integer v [for a momentum independent irreducible part or, equivalently, 
self-energy this statement follows directly from (jH]) and ([7])] . In view of this re¬ 
sult it is convenient to split the first Brillouin zone into n streaks of the width 
which are oriented parallel to the x axis. Let us denote wave vectors in the 
streak with —-<qy<—- + — and — - < Oa, < - as k. Then any wave vector 
in the first Brillouin zone can be represented as h + jK with 0 < j < n — 1. In 
these notations 


G'’"(q'qm) = G'^(k + j'K,'k + jK,,m) = GJ,j(km). 
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Analogous notations can be used for K'^{q'qm) and r‘^(q'q). Considering in¬ 
dices j and j' as matrix indices equation ([5]) can be rewritten as 

O'"(km) = [1 - K'"(km)T'"(k)]”^ K'"(km), (8) 

where matrices are denoted by the boldface font and 1 is a n x n identity matrix. 

If the approximation of the dynamic mean field theory |82] - a momentum 
independent self-energy - is accepted, Eq. ([8]) acquires the form 

G-(kcc) = {[if'^(a;)]-'l-T'^(k)}”\ (9) 

where the analytic continuation to the real frequency axis was performed. In 

this approximation each Landau subband forms strongly correlated bands in¬ 
dependently of other subbands. Indeed, let us denote the part of T'^(k), Eq. 

which is proportional to t, as T'(k). It is a Hermitian matrix with the 
eigenvectors VA(k) and eigenvalues ^^(k), 

T'(k)VA(k) =EA(k)VA(k), 0<A<n-I. (10) 

This equation is the Harper equation [33] for calculating dispersions ^^(k) of 
the Landau subbands. The vectors VA(k) are also the eigenvectors of the matrix 
G'^(kw), 

G"(kLc)VA(k) = 

xVA(k) (II) 

Except the Zeeman contribution the braces in the above equation look like the 
poles of the electron Green’s function of the Hubbard model with the initial 
dispersion E\(k) in zero field [23 [52]. Hence each Landau subband A forms 
strongly correlated bands independently of other subbands. In other words, 
with respect to strong correlations the Landau subbands behave as independent 
bands in this approximation. Notice also that for vanishing U only the first- 
order cumulant remains nonzero, [K'^{uj)] ^ ^ oj — fj, and the above equations 
reduce to Green’s function of uncorrelated electrons in the magnetic field. 

From Eq. (fTTll we find for the DOS 

p'^iu;) = -^Im^TrG'^(kc.) 

k 

kA ^ 

1 1 “' 

, (12) 

where N = nN^Ny is the number of sites. 

In the below calculations the irreducible part K'^{ui) was approximated by 
the first-order cumulant C'i(t't) = (TdiCT(T')aiCT(T))o, where the subscript 0 
indicates that the averaging and time dependencies are determined by Hg. In 
the absence of the magnetic field this approximation leads to the Hubbard-I 
approximation [21j . For the chemical potential satisfying the conditions 

T</i, T-^U-p, (13) 
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the cumulant reads 


(14) 


Ciiio) = 


LU IJ, -2 


(_UJ fl — U) 

Substituting this approximation for the irreducible part into (fT^ we find 

,'1 _ _ I _ 

^ ~ ‘2N ^ v^C/2 +£'2(kACT) 

X { \yTpTWi^^ + ^(kAa)] <5 (w - 4;,+) 

+ [x/[/2+£:2(kAa) - ^(kAa) 


(15) 


where 

^(kAcr) = Exik) + 

(16) 

4a± = ^ + £’(kAcr)] ± + £^{k\a) - n 

and Ex (k) is a solution of the Harper equation (ITUl) with the tridiagonal cyclic 
matrix T'(k) from ([T]). 


3 Results and discussion 


Below we consider the density of electron states (USD. For T = 0 the integral 
of the product of this quantity and the frequency over the occupied states gives 
the electron contribution to the thermodynamic potential H. As in the case of 
weak electron correlations, oscillations of p{uj = 0), the DOS at the Fermi level, 
in varying magnetic field lead to oscillations in H and its derivatives, which are 
observed in quantum oscillation measurements m- 

Let us first consider the difference in the oscillations of p'^(O) for electrons, 
which spins are parallel and antiparallel to the applied magnetic field. The 
difference is caused by the Zeeman term of Hamiltonian (ED- An example is 
given in Fig. [T] In these calculations we have substituted (5-functions in m 
with Lorentzians with the artificial broadening 77 = 0.03t, which imitates finite 
lifetimes of states. The oscillations are very sensitive to this broadening: with 
increasing 77 the oscillation amplitude decreases rapidly. The parameter z = 
indicated in the figure caption, characterizes the value of the Zeeman 
term. For a = dA, the intersite distance approximately corresponding to the 
distance between Cu ions in Cu-0 planes of YBa 2 Cu 307 _y, this parameter is 
equal to 1.5 eV. For zero temperature the electron concentration is calculated 
using the formula 


n = 





(17) 


As seen from the figure, the oscillations in p*^“^)(0) and p*^+^^(0) are shifted in 
phase, and this shift grows with growing zjt. However, the frequencies of the two 
oscillations are the same and equal to F 1.7 kT. Notice that this frequency by 
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Figure 1: (Color online) The densities of states at the Fermi level (the blue 
solid line), (the red dashed line) and their mean value ^ 

(the olive dash-dotted line) as functions of the inverse magnetic induction ex¬ 
pressed in terms of ^ T = 0, U = 8t, n' — 3, fi = 0.953 and = 2<. 


an order of magnitude smaller than oscillation frequencies in metals with large 
Fermi surfaces |13j . The frequency does not change with the variation of B or 
z. Due to the phase shift the amplitude of oscillations in p = ^ + P^'*’^^)) 

which enters into D, is smaller than in any of the two summands. However, 
even in the worst case when oscillations in them are in antiphase, as in Fig. [TJ 
the complete compensation of oscillations does not occur in p. 

As follows from the symmetry of the Hamiltonian, for w = 0 

p‘"(t, p) = p“'"(t, H - p) = p'^(-t, p). 

Therefore, below we set t > 0 and p < ^, which corresponds to h < 1. Be¬ 
sides, we neglect the Zeeman term, which does not qualitatively change results 
discussed below. The variation of the DOS oscillation with the electron filling 
is shown in Fig. [2l In the figure caption, the oscillation frequencies F are esti¬ 
mated by setting a = dA, as above. As seen from the figure, the frequency F 
is changed by an order of magnitude when the electron filling is decreased from 
n = 0.969 to 0.751. This corresponds to the shift of the chemical potential from 
the value p = 1.4t near the top of the lower Hubbard subband to p = 0.4t. 

What is the reason for small values of F at small deviations from half-filling? 
As was indicated in the previous section, a Landau subband forms strongly 
correlated bands independently of other subbands. From (ITCl) we find for the 
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Figure 2: (Color online) The density of states at the Fermi level as a function of 
the inverse magnetic induction expressed in terms of ^ ^ “ 0.969 

(a, F = 0.95 kT), n = 0.897 (b, F = 3.9 kT) and n = 0.751 (c, F = 8.6 kT). 
T = 0, U = 8t and n' = 3. 


energies of Landau subbands contributing to the DOS at the Fermi level 

i5.(k) . (18) 

This equation is obtained for the case h < 1 ignoring the Zeeman contribution. 
Notice that for /r —>■ 0, when the chemical potential goes deep into the lower 
Hubbard subband and the value 1 — h is large, i?A(k;) —> 0. Landau subbands 
ExCk) obtained from the Harper equation (1101) are located in the energy range 
(—Y, ■§■), where A = 8< is the width of the initial band. The above result means 
that for large deviations from half-filling Landau subbands near the center of this 
range contribute to /o(0). In this spectral region the number of Landau subbands 
is large and their positions change rapidly with the variation of B dmiH]. This 
leads to high-frequency oscillations of /^(O) with changing induction for such 
electron concentrations. On the other hand, the value of the chemical potential 
near the top of the lower Hubbard subband, when 1 — h is small, can be also 
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estimated from (USD, 



Substituting this result into (fT51) we get i?A(k) ~ other words, in this case 

Landau subbands at the periphery of their energy range contribute to p(0). The 
number of these subbands is smaller and the distance between them is larger 
than in the central part of the Landau spectrum [221 [28]. This results in the 
slower frequency of oscillations observed in p(0) for these 

The above discussion is valid for the case U > A. Though the Hubbard-I 
approximation used above is rather rough for smaller U, we hope it gives at least 
qualitatively correct results. For t7 <C A widths of the bands ekA±, Eq. (117)1) . 
become larger, and a larger number of the Landau subbands contributes to p(0). 
As a result, in contrast to the case of strong correlations, for /r ~ Utop the main 
contribution to p(0) is made by subbands, which are well off the spectrum edges 
±A. Therefore, for [/ <g; A the variation of the oscillation frequency with the 
change of the electron filling is less pronounced in comparison with the case 
[/ > A, and the oscillations are observable only for large field intensities. 

For U > A the frequencies F obtained above for small deviations from 
half-filling are of the same order of magnitude as those found in lightly doped 
cuprates [i nni HB (nj. The mechanism, which leads to these decreased fre¬ 
quencies, is based on the important feature of these crystals - strong electron 
correlations. This gives promise that the mechanism can be used for the inter¬ 
pretation of the mentioned experiments on quantum oscillations. 


4 Concluding remarks 

In this work the two-dimensional fermionic Hubbard model in a perpendicular 
homogeneous magnetic field was considered. Using the strong coupling diagram 
technique expressions for the electron Green’s function and the density of states 
were derived. Calculations were performed for the nearest neighbor form of 
the kinetic energy and for the approximation of the irreducible part, which 
corresponds to the Hubbard-I approximation. For a momentum independent 
irreducible part/self-energy each Landau subband was found to form strongly 
correlated bands independently of other subbands. The density of states at 
the Fermi level as a function of the inverse magnetic induction demonstrates 
oscillations, which lead to oscillations in the thermodynamic potential and its 
derivatives, observed in quantum oscillation measurements. Due to the Zeeman 
term the densities of states of the electrons with spins parallel and antiparallel 
to the magnetic field are shifted in phase. However, the full compensation of 
these oscillations in the combined density of states does not occur. For Hubbard 
repulsions comparable to the width of the initial electron band frequencies of 
oscillations for small deviations from half-filling are much smaller than those 
for larger deviations. The reason of this is the fact that for different electron 
fillings different parts of the spectrum of Landau subbands contribute to the 
density of states at the Fermi level. For small deviations from half-filling it is 
the periphery of this spectrum, while for larger deviations it is a part near the 
center of the spectrum, which varies faster in changing field than the periphery. 
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Oscillation frequencies for small deviations from half-filling are comparable with 
those observed in quantum oscillation experiments on lightly doped yttrium 
cuprates. 

This work was supported by the research project IUT2-27, the European Re¬ 
gional Development Fund TK114 and the Estonian Scientific Foundation (grant 
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